Mutual coherence of optical and matter waves 
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We propose a scheme to measure the cross-correlations and mutual coherence of optical and 
matter fields. It relies on the combination of a matter-wave detector operating by photoionization of 
the atoms and a traditional absorption photodetector. We show that the double-detection signal is 
sensitive to cross-correlation functions of light and matter waves. 
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I. INTRODUCTION 



II. THE DETECTION SCHEME 



The manipulation and control of matter waves by op- 
tical fields forms a new paradigm that is now receiving 
much attention. In a first generation of experiments, ex- 
ternally applied optical fields were used to manipulate ul- 
tracold atomic samples. Examples include optical dipole 
traps for multicomponent condensates B and the diffrac- 
tion of condensates by optical gratings ||,|| . More recent 
experiments are also concerned with the back-action of 
the condensate on the light field. For example, Bragg 
scattering of light off a condensate has been used to de- 
termine its form factor and recently, superradiance 
experiments [gj have demonstrated the joint build-up of 
optical and matter- wave fields. One topic of particular 
interest in this context is the possibility to generate cross- 
correlations and quantum entanglement between optical 
and matter waves. Indeed, there are already several the- 
oretical predictions of such effects [|]-|9|, e.g. in an ultra- 
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cold version of the Collective Atom Recoil Laser [[Hi 

It is therefore timely to investigate ways to measure 
the joint quantum statistical properties of optical and 
matter- wave fields. This is an extension from the usual 
quantum optics situation, where one normally measures 
either the statistical properties of light, or more rarely, 
those of the atoms, but not joined properties of the cou- 
pled system. For this latter purpose, it is necessary to 
develop a generalized detection theory applicable to the 
simultaneous detection of optical and matter-wave fields. 

Our approach builds on the well-known results of opti- 
cal detection and coherence theory ||l2| and their recent 
extension to atomic matter waves |13| . Specifically, we 
consider a detection scheme consisting of an atomic de- 
tector operating by ionization [T^ ] and a traditional pho- 
todetector operating by absorption ||T^ , |l4j| . Our model 
is introduced in section II. Section III presents general 
expressions for the correlations between optical and mat- 
ter fields. Our general results are illustrated in section 
IV, where mutual coherence functions are calculated for 
a simple model of linear coupling between optical and 
matter waves. Finally, section V is a summary and con- 
clusion. 



We consider a situation where an optical field and an 
atomic matter-wave field are dynamically coupled, and 
we want to quantify the cross-correlation and mutual co- 
herence of these fields. The matter- wave field might, but 
needs not to be, a Bose-Einstein condensate, described 
by a multicomponent field operator ^(r), satisfying the 
bosonic commutation relations 



[*i(r),V(r')] 



5ij8(r - r'), 



(1) 



where r is the center-of-mass coordinate of the atoms 
and the indices i label their internal state. The optical 
field is described as usual by the field operator E(r, t). 
The dynamics of the combined system is governed by 
the Hamiltonian 



Ho = TL S + H. r , 



V 



(2) 



where 7i s and TL m describe the evolution of the 
Schrodinger atomic matter- wave field and of the Maxwell 
light field, respectively, while V is the interaction between 
these two fields, typically the electric dipole interaction. 
It is not necessary to give more details about the specific 
system at hand at this point. 

We now turn to the description of the detector. It con- 
sists of two parts: The first one is a conventional optical 
photodetector that emits photoelectrons when absorbing 
photons from the light field. In addition, there is also 
a matter- wave detector. It consists of a tightly focussed 
laser beam that can ionize atoms, thereby annihilating 
them and producing a ion-electron pair. The detection of 
the mutual coherence properties of the coupled Maxwell- 
Schrodinger fields is obtained from the correlations be- 
tween the electron emitted by these two detectors. 

Following the original work of Glauber jL2j, we de- 
scribe the photodetector as consisting of two-level atoms 
whose excited state is in the continuum, corresponding 
to photoionization. We call it the Maxwell detector. Its 
interaction with the system is decribed as usual by the 
electric dipole interaction Hamiltonian 



V m = -d • E(r m ), 



(3) 



where d is the atomic electric dipole moment between 
the detector ground electronic state and the continuum 
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state, and r m is the location of the detector, assumed to 
be punctual for simplicity. 

The interaction between the system and the 
Schrddinger detector, also assumed to be punctual at lo- 
cation r,, is 



V s = ?i^^(r s )[*i(r s )*(r s )e-^* + h.c.\. 



(4) 



Here ^> and 'P, are the field operators associated with 
atoms in the sample and the continuum state \i), respec- 
tively, tti is a coupling constant proportional to the elec- 
tromagnetic field producing the ionization and lol is the 
laser field frequency. The assumption of a punctual pho- 
todetector is well justified for ultracold atomic samples, 
in which case the ionizing laser can easily be focussed 
to a spot much smaller than a characteristic atomic de 
Broglie wavelength. 

The Hamiltonian describing the total interaction be- 
tween the detector and the system is 



v d = v m + v s 



(5) 



and the mutual coherence of the optical and atomic fields 
is determined to lowest order by measuring the joint 
counts of the Maxwell and Schrodinger detectors, i.e., 
the joint probability to count one electron 



III. CORRELATIONS BETWEEN DETECTED 
SIGNALS 

To lowest order, the probability amplitude for the joint 
detection of electons from the Schrodinger and Maxwell 
detectors is given by second-order perturbation theory 
in Vd- To that order, the transition probability from the 
system initial state \i) to a final state |/} in a time interval 
At is 



1 



t+At 



dt" 



t+At 



dt'(f\Vd(t")V d (t')\i), (6) 



where t" > t', and the detector interaction Vd is now 
expressed in the interaction picture 



Vd(t) 



J(H +H d )t 



V d (t)e 



-i(H a +n d )t 



(J) 



TLd being the detector Hamiltonian. 

When it comes to measuring the cross-correlations be- 
tween the optical and atomic fields, the relevant terms are 
of course those involving one electron each emitted by the 
Schrodinger and Maxwell detectors. In Eq. (JsJ) , these are 
the terms involving cross-products of the detector-system 
coupling Hamiltonians V s and V m , that is, 



\ / dt" f dt' 

h Jo Jo 

[(f\V m (t")Vs(t')\i) 



(f\V s (t")V m (t')\i)}, (8) 



with t" > t'. 

The probability to jointly excite the Maxwell and 
Schrodinger detectors is therefore 

P(Ai)~-^Vft m (/ m )7^(/ s ) 



t+At 



dt" 



t+At 



dt' {(f\V m (t")V s (t')\i) 



{f\V s (t")V m (t')\i)}\ 2 , t">t'. 



(9) 



The sum is on a complete set of final states {/} = {|/o) <8> 
|/m) <$ \fs)}, where {|/o}} is a complete set of states of 
the system, and lZ m (f m ) and lZ s (f s ) are the sensitivities 
of the Maxwell and Schrodinger detectors, respectively, 
associated with the final states |/ m ) and |/ s ) of these 
detectors. This sum results from the fact that we are not 
interested in the final state of the system. It allows one 
to introduce a closure relation and to reexpress Eq. (|^) 
as 



^ rt + At rt + At rt + At p 

P(At) ~ -2 / dh dti / dt 3 \ 

n Jt Jt Jt Jt 

x [rj s (r s ,t4,t 3 )T] m (r m ,t2,ti) 

x (E- (r m , t 2 )¥ (r s , t 4 )tf (r a , t 3 )E+(r m , h)) 

+i] s (r s ,t 1 ,t 2 )T] m (r m ,t4 : ,t 3 ) 

x (&(r s ,t 2 )E- (r m ,U)E+{r m , t 3 )tf (r a , h)}] 



rt+At 



(10) 



where 



Vs(r s ,t,t') = ie^ ]T]T 7e s (/ s )|^(r s 

i fs 

x(*.|* i (r a ,t)|/.>(/.|* i t (r a ,OI«»> 



(11) 



and 
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r/ m (r m ,t,t') = ~^U m (f m )(i m \d(t)\f m )(f m \d(t')\i m ) 

fm 

(12) 

are the detector sensitivities, and d(t) is the atomic dipole 
operator in the interaction picture. 

We have neglected contributions involving the product 
of two creation and annihilation operators. This follows 
from the implicit assumptions that the Maxwell detector 
is initially in its ground state, and that in the Schrodinger 
field detector any atom in the continuum will be removed 
from the sample instantaneously. Note also that the sta- 
tistical properties of the atoms in the detectors do not 
come into play jljj, a consequence of the fact that the 
lowest order cross-correlations that we consider here in- 
volves only one measurement for each of the fields in- 
volved. 

Equation ( |l0| ) can be significantly simplified in the case 
of broadband detection, that is, for situations where the 
bandwidths of the detectors are much broader than those 
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of the detected fields. In this case, and for stationary pro- 
cesses, the Markov approximation can be invoked in the 
treatment of the detectors. Their responsivities can then 
be approximated as <5-functions of the time differences 
t% — t\ and t± — <3. 

Before proceeding with a final expression for the joint 
counting rates, we observe that expression (10) involves 
the sum of two normally-ordered correlation functions. 
This is a consequence of time ordering and of the non- 
commutivity of the Schrodinger and Maxwell field oper- 
ators at different times. One possible way to measure 
each of these correlation functions separately is by using 
a time-gated detection scheme |14 with nonoverlapping 
in time windows for the two detectors. Specifically, if the 
ionizing laser beam of the Schrodinger detector is turned 
off after a time t s and the photodetector is masked af- 
ter t m , then the coupling of the detectors to the system 
becomes 



V = V m 0(t m -t)+V s 9(t s -t), 



(13) 



where we can chose for example t s < t m < At. In this 
particular gating scenario, only the first term in Eq. (|Io| ) 
remains. 

Assuming that the spectrum of the Schrodinger and 
Maxwell fields to be analyzed are centered at ui a and u) m , 
respectively, the probability of registering one count at 
each detector during the time intervals t s and t m is then 

Wm)f?s( r s, U s ) / dt 2 / dti 

Jo Jo 

x (E^(v m , t 2 )*Hr s ,t 1 )V(r s ,t 1 )EW ( Fm) h )) (14) 

where i] m (r m ,Lu m ) and rj s (r s ,oj s ) are the Fourier trans- 
forms of the detector Maxwell and Schrodinger respon- 
sivities. Hence the counting rate is simply propotional to 
the joint correlation function of the fields, 



d 2 



dt s dt, 



-V(t s ,t m ) = r) m (r m ,Lj m )r} s (r s ,u s ) 



x (E ( -\v m ,t m )^(r s ,t s )^{v s ,t s )E^(r m ,t m )). (15) 

For t m < t s < At, it is only the second term in Eq. 
( |i"o| ) that persists. It is straightforward to modify the 
final form ( [15|) o f the counting rate accordingly. 

Equation Q15|) is similar to the double-photocount re- 
sult in conventional optics. This is of course not an ac- 
cident, since the only difference is that instead of hav- 
ing two photodetectors, we have replaced one of them 
by a matter-wave detector. However, there is an es- 
sential difference since this new detection scheme allows 
one to measure joint coherence properties of the optical 
and matter fields. In addition, it is obviuous that the 
proposed detection scheme can readily be generalized to 
measure higher-order correlatinos of these fields. In that 
case, though, the fermionic quantum statistics of the ions 
and electrons resulting from the photoionization of the 
bosonic matter- wave field must be handled properly [|l3| . 
This will be the subject of a future publication. 



To see how the lowest-order cross-coherence measure- 
ment of Eq. ([l^) works in practice, we illustrate it for 
the case of a simple example. 



IV. EXAMPLE 

The simplest situation one can consider is that of a lin- 
ear coupling between the optical and matter-wave field. 
Of course, we should keep in mind such a coupling is not 
possible in principle, due to the conservation of the num- 
ber of particles. Instead, all couplings must be at least bi- 
linear in the Schrodinger field operators, so that an atom 
annihilated in one state is then created in another state. 
However, the advent of Bose-Einstein condensates with a 
large number of atoms in a single, macroscopically popu- 
lated quantum state, allows one to describe that state as 
a c- number. This is analogous to the undepleted, classi- 
cal field approximation in quantum optics. In that case, 
the bilinear matter- wave part of the interaction Hamilto- 
nian, which involves both the condensate mode and some 
other side mode, approximately reduces to a linear form 
only (see e.g. P,pj[). With this kind of situation in mind, 
we then consider a matter-light interaction Hamiltonian 
of the general form 



N 



M 



H s = fttJ a a) a a a + hujiclci 

a—l i—1 

+ [5i 1)a a Ci + 9i 2) a a Ci + H.c. 



(16) 



where N and M are the number of modes in the optical 
and matter field, respectively, with eigenfrequencies w a 
and uii. (We use greek indices for the optical field modes 
and roman letters for the matter waves in this example.) 
The mode annihilation operators 



and 



dru*(r)E+(r) 



drcf)*{v)^(v) 



(17) 



(18) 



are bosonic optical field mode operators, u Q (r) and 4>i(r) 
being the corresponding mode functions, E + (r) is the 
positive frequency part of the electric field operator, 
9i ■> 9^ are coupling constants, which allow for a para- 
metric amplification type of matter-light coupling as well. 
These constants are proportional to the order parameter 
of the condensate ground state, as discussed above. 

In addition, we allow for both the optical and the mat- 
ter modes to be coupled to thermal reservoirs which re- 
sult in losses at rates K a and Ki. The Langevin equations 
describing the evolution of this model system are n& 



-x(t) = -Mx(t) + B£(i), 



(19) 



3 



where x(t) = {a a (t), a(t)} T , a = {l,...,N};i = 
{1, . . . , M}, and M is a (N + M) x (N+M) matrix whose 
elements depend on the system parameters. The diffu- 
sion matrix in Eq. (|l9| ) is D = BB^ = diag(n a K a ,riiKi), 
where n a and fii are the thermal populations of the reser- 
voir modes. We assume for simplicity that the reservoir 
operators can be treated in the Markov approximation 
as 5-correlated white noise with zero expectation value, 
i.e. (&(i)> = and (&(t)fc(f )) = M(* ~ *')-. 



For this example, the counting rate of Eq. ( 15 ) reduces 
to a sum over all possible fourth-order cross-correlation 
functions 

W(r s t s ,r m t m ) = ^ u a{ V rn)Ui3{Y m )4> k i {Y s )(j) J {Y s ) 

a/3ij 

+ (U-^G(t s ,t s )(U-^).. (U- l G{t mi t m ){U- i y) ap 



where U is the matrix of eigenvectors of the system ma- 
trix M, so that 



U~ 1 M.U = diag{\ k ), k = {l,...,N+ M}, 
where {Ai} are the eigenvalues of M, and 

e -Xi(t-s) 



(21) 



G tJ (t,s) = (C7D[/t) 
Gij (t, s 



11 + 

p -A*( S -t) 

(UDU% 



A i + A* 



t > s, 



t < s. (22) 



The two terms in square brackets in Eq. ( |20|) correspond 
to the stationary value of the counting rate and to the 
exponential approach to that value. These expressions 
follow directly from the stationary solution of Eq. (|l^) , 



x(f) 



,-M(t-t'] 



B£(t')dt', 



(23) 



where x(0) = {a a ,Ci} T . 

The possibility to extract particular mode correlation 
functions from these general expressions depends on the 
specifics of the problem. To gain some insight, we con- 
sider the two-mode problem with the drift matrix 



M 



and diffusion matrix 



D 



Ki/2 

ig 











(24) 



(25) 



g (2) (T) 



1.0 




0.0 

0.5 1.0 1-5 ^ 

FIG. 1. Time dependence of the normalized station- 
ary cross-correlation function g, 
ter-wave modes for g = 10., Ki 
n 2 = 0.1 



^ (r) of the optical and mat- 
: 6., K2 = 6., ni = 0.1 and 



(20) 



The normalized stationary cross-correlation function 
gx (t) of the optical and matter-wave modes is shown 
in Fig.l as a function of the time difference t = t m — t s . 
As expected, for this case of linear system with posi- 

(2) 

tive diffusion coefficients, g x (t) exhibits bunching and 

(2) 

lim T _>cx3 gx (t) = 1. The oscillations correspond to the 
fact that for the chosen parameters the eigenfrequencies 
of the system are complex, with a real part smaller than 
the imaginary part. 



V. DISCUSSION AND CONCLUSIONS 

In this paper, we have generalized the familiar ap- 
proach to optical detection and coherence theory to 
deal with situations involving coupled Maxwell and 
Schrodinger fields. We have proposed a concrete detec- 
tion scheme that allows one to experimentally determine 
cross-correlation functions and the mutual coherence of 
these fields. We have concentrated only on the lowest- 
order cross-correlation, but the scheme is obviously eas- 
ily generalized to measure higher-order correlations. In 
that latter case, though, proper care must be taken of 
the quantum statistics of the detected signal. This topic, 
where the fact that atoms are actually composite bosons 
(or fermions) plays a central role, will be discussed in a 
future publication. 

We have illustrated our detection scheme on a simple 
case of linear coupling between an atomic matter wave 
and an optical field. The results in that case are quite 
similar to the situation involving only optical fields. This 
is a direct consequence of the linearization procedure 
leading to our model Hamiltonian. More realistic situ- 
ations will require the use of an interaction which is at 
least bilinear in the matter-wave field, in which case ma- 
jor differences can be expected, for instance in situations 
involving small condensates. 

The detection scheme presented here, as well as its 
generalizations, are expected to find applications in a 
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number of situations involving e.g. the optical control 
and manipulation of matter waves, nonlinear atom op- 
tics, the generation and detection of quantum entangle- 
ment between light and matter, and other applications 
of quantum atom optics. Note however that some of our 
simplifying assumptions, such as e.g. the model a punc- 
tual detector, only hold provided that the ionizing laser 
can be focused to a spot small compared to an atomic 
de Broglie wavelength. As such, it is geared toward ul- 
tracold samples. It would not be difficult to expand our 
scheme to detectors of finite size, but in that case, the 
detected photocounts will involve a convolution with the 
spatial detector resolution, with a concomitant washing 
out of the finer spatial structure of the cross-correlations. 
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